CCNY-HEP-06/7 
May 2006 



Color Skyrmions in the Quark-Gluon Plasma 



O 

^ ; JIAN DAIi and V. P. NAIR^ 

, Physics Department 

City College of the CUNY 
O ; New York, NY 10031 

(N 
> 
O 
0^ 
O 

g ■ Abstract 
O , 

I We consider the general formulation of nonabelian fluid dynamics based on symmetry con- 

QhI siderations. We point out that, quite generally, this admits solitonic excitations which are 

Qh! the color analog of skyrmions. Some general properties of the solitons are discussed. 

> 

X 



^E-mail: jdai@sci.ccny.cuny.edu 
^E-mail: vpn@sci.ccny.cuny.edu 



1 Introduction 



The collision of heavy nuclei at the Relativistic Heavy Ion Collider has created an interesting 
new state of matter [1]. It is clear that a state with deconfinement of quarks and gluons 
has been achieved. There are various indications that the resulting fluid is probably best 
described as a color liquid. One of the surprises has been the very low shear viscosity of the 
fluid. The value of the viscosity, it seems, is close to what may be the theoretical lower limit 
possible, a value which may be understandable in terms of a gravity-dual to the theory [2]. 
Thus a good first approximation to the description of this color liquid may be as a 'perfect 
fluid'. However, unlike the case of ordinary fluids, since the constituents carry color degrees 
of freedom, the transport of such degrees of freedom, in a way consistent with the nonabelian 
gauge symmetry, becomes an important issue for this 'color liquid'. This can, of course, be 
studied via kinetic equations, starting from the basics of QCD. However, the hierarchy of 
kinetic equations has to be truncated for reasons of computability, very often to the level of 
uncor related single-particle distribution functions, and such an approach is then limited to 
dilute systems near equilibrium. Since the experimental results indicate that the color liquid 
is not a dilute system, the validity and fruitfulncss of this approach become questionable. By 
contrast, experience with ordinary fluid dynamics shows that one can derive the equations 
of fluid dynamics from very general principles, which then shows that the equations have a 
range of validity significantly beyond the regime where the truncated kinetic equations apply. 
This latter, a priori, approach has been developed for nonabelian fluid dynamics as well [3, 4]. 
The nonabelian degrees of freedom of the fluid are described by a field which takes values in 
the color group, for example, by g{x,t) G SU{3) for QCD. An immediate and qualitatively 
striking consequence is that nonabelian fluid dynamics leads to topological solitons. These 
are color skyrmions associated with the homotopy group n3[6'C/(3)] = Z. The topological 
quantum number is the color analog of what has been known as helicity in ordinary fluid 
dynamics for a long time. The purpose of this article is to point out the existence of such 
solitons, or configurations of nonzero color helicity, and study some of their properties. 

In the next section, we will review briefly the formulation of ordinary fluid mechanics in 
group theory language and its generalization to include the transport of nonabelian degrees 
of freedom. This subject has been reviewed recently in [5]. In section 3, we will introduce 
the solitons and work out their properties. We then conclude with a short discussion. 

2 Group theory and fluid dynamics 

We begin with the well-known observation that ordinary fluid dynamics can be described 
as a Poisson bracket system. With p as the fluid density and Vi as the fluid velocity, the 
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fundamental Poisson bracket relations are given by 



[ p{x,t),p{y,t)] 
[vi{x, t),p{y,t)] 



0, 



(1) 



where Wj,- is the vorticity defined by 



ujij = diVj - djVi 



(2) 



Equations (1) will lead to the usual equations of fluid dynamics, where the Hamiltonian H 
may be taken as 



d-^x 



^pv^ + V{p) 



(3) 



The canonical structure, or the symplectic form, defined by the Lagrangian is the inverse of 
the fundamental Poisson bracket. In the present case, the fundamental Poisson brackets have 
a zero mode and so this is a complication in finding a suitable Lagrangian. This can be seen 
in terms of the Chern-Simons action for the velocity, namely, 



8ttJ 



VidjVk 



(4) 



(The invariant defined by C is the fluid helicity.) It is easily seen that C Poisson commutes 
with any observable F, 



5F_\ 

_ Sp dx'^ \ Svi J Sp dx^ \ Svi J 



SF d f 5C \ _6C d 



- - W,; 



S£SC_ 
' Svi Svj 







(5) 



This shows that SC/Svi is a zero mode for the brackets, preventing the inversion needed to 
deflne the action. The solution to this problem is also rather clear. Since C commutes with all 
observables, we must fix its value and consider only variations in Vi which preserve this value. 
In particular, if we choose C = 0, we must consider Vi which trivializes the Chern-Simons 
term by making VidjV^e^^'^ into a total derivative. Such Vi are given in terms of the Clebsch 
parametrization, 

Vi = die + a diP, (6) 

where 9, a, f3 are arbitrary functions. 

A suitable action which leads to the Poisson bracket relations (1), and the continuity and 
Euler equations, is 



S = j dtd^x 

Notice that we have the canonically conjugate pairs (p, 6), {pa,P) 



pe + pap-^pv^- V{p) 
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We will now rewrite this in terms of a group element 

^ " x/r^ \ u 1 ; \ ; 

where u is a complex variable, g is an element of SU{1, 1) if u = u and an element of SU{2) 
if {t = —u. (We will consider both possibilities together for a while.) Using (8), we find by 
direct computation that [6, 3] 

-iTr{a3g~'^dg) = 6.6 + adp (9) 

where 

a = ——: , p = --log{u/u) (10) 

J- Utl Zi 

Equations (9,10) show that the Clebsch parametrization can be expressed in terms of a single 
group element as Vi = —W^{o-ig~^ d%g\ Notice that Q in (8) corresponds to the (T3-direction 
in g\ thus a, /3 parametrize the space 5*^7(1, 1)/C/(1) and S'?7(2)/(l), respectively, for u = u 
and u = —u. Going back to (7) we notice that the action can be written as 



(11) 



where j° = p. Elimination of f using its equation of motion takes us back to (7). The form 
of equation (11) also shows us how to generalize to a relativistic situation, 

S = j d'^x [-if Tiia^g-^d^g) - F(n)] (12) 

where = ffg^v, g^v being the metric tensor. The choice of the function F specifies the 
equation of state for the fluid. A four-velocity u^^ for the fluid can be defined by f = nu'^. 
The energy-momentum tensor is then 

T"" = nF'uf'u'' - g'^'inF' - F) (13) 

which also identifies the pressure as p = nF' — F. 

The Poisson bracket relations which follow from (11) (or (12)) are 

[p{x,t),p{y,t)] = 

[p{x,t), g{y,t)] = -ig{x) ts S^^\x - y) (14) 

Here ta = ^aa- This equation shows that p generates right translations on g along the 
(T3-direction. In the quantum theory, it is easily checked that U = exp[— 27ri / p] acts on g as 

U^gU = g e^""'*^ = -g (15) 
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All observables are invariant under U since they have even powers of g, This means that 
U is the identity operator on observables, which, in turn, implies that J p = integer [7]. This 



is essentially the statement that the fluid is made of particles. Mathematically, this is due 
to the ^-direction being compact. Turning this logic around, we see that it is consistent to 
require that the field 9, which appears in the Clebsch parametrization, should be a compact 
direction since fluids of interest are ultimately made of particles. 

Prom the Clebsch parametrization, coij = (diadj P — djadiP) . Consider now the integral of 
this over some compact region V in M^. Jyoo is the volume of 5*^7(1, 1)/C/(1) or SU{2)/U{1) 
over the image of V via the map g{x). This will be quantized for S*^ = SU{2)/U{1), but not 
for SU (1, 1)/?7(1). Thus we see that the choice of the group is determined by the quantization 
of J p and by whether we want quantized vorticity or not. 

We now turn to the generalization of this to other groups. A Clebsch-type parametrization 
for higher groups would be an obvious direction to try [3]. This is possible for the groups 
SO{n) with = —iTi{t°'R^diR) where R G S0{2n — 1) and ta are generators of SO{n) C 
SO{2n — 1). The relevant coset is S0{2n — 1)/S0{n) x SO{n — 1). The resulting dynamics 
is rather involved and does not admit the Eckart factorization = Q°'ji which is what we 
would expect if the nonabelian charge density Q"- is transported by particle motion. 

A different approach is to start with the motion of particles carrying nonabelian charges. 
This is described by the Wong equations [8] 



where Q"" = TT{a^g~^t°'g) and g G SU (2). (Here we are considering SU{2) color for simplicity. 
We will extend this to any Lie group shortly.) The first term is the usvial action for a 
particle coupled to an external gauge field except for the color charge factor Q". The second 
term is what leads to the dynamics of the color charge. Notice that is invariant under 
g — >■ gh, where h = exp(it3(^). The second term in the action changes by a surface term 
AS = n / ifdt = nt^ip. Thus, the equations of motion are invariant under g gh. However, 
one can consider closed loops in the [/(l)-subgroup, generated as a trajectory over time, 
which would have A(p = 2it. The invariance of e*'^ then requires n E Z. For integer values 
of n, the theory (including the quantum theory) involves only SU{2)/U{1) = S'^ variables. 
In fact, the second term of S gives the symplectic form = inTr{a3g~^dg A g~^dg) for 
the parameters in g. Since the phase space for this is the two-sphere, the phase volume is 



qa _ jabc^c Qb ^ Q 

Pi-Ft.x^Q'' = 



(16) 



These equations of motion can be obtained from the action [9] 




(17) 
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finite, giving a finite-dimensional Hilbert space upon quantization. In fact, the quantization 
of (17) leads to one unitary irreducible representation of SU{2) with a highest weight state 
of j = n/2. The action (17) thus leads to the standard description of color by matrices. The 
term —inTr{a3g~^dg) is often referred to as the Kostant-Kirillov-Souriau (KKS) form [10]. 

Focusing on the color degrees of freedom, the generalization of (the second term of) the 
action (17) to many particles is given by 

S = -i fdtj^nx 1V(a3 g^^gx) (18) 

where A = 1,2,. . . ,N, labels the particle under consideration. We can now take the contin- 
uum limit by taking N large, with an average volume v per particle. This means that we can 
write X ^ X, J2^^ J d^x/v and nx/v = Jq. The action (18) then becomes 



I J d^x f T^ia^g-^g) (19) 



where g = g{x,t). Notice the similarity of this to the first term of the action for ordinary 
fluid dynamics. Taking (19) as the key term which leads to the symplectic form for the color 
degrees of freedom, we can write the action for fluids, where the particle carry SU (2) color 
degrees of freedom, as 

1 



d'^x 



-if Tr{a3g-'D^g) - F{n) - -F^.F-'^^ 



(20) 



where D^g = d^g + A^g and = —it"'A^^ is a nonabelian background field. (By this we 
mean the field background in which the fluid moves; thus A^ represents gluon degrees of 
freedom which are not homogenized with the fluid, such as hard gluons.) j'* is a current 
four- vector, = j'^fguv f can be eliminated by its equation of motion to get a simpler 
form of S, which, however, will not be manifestly Lorentz invariant. The dynamics which 
follows from (20) will be the SU (2) analog of magnetohydrodynamics. The color current, 
which couples to A^^, is easily seen to be 

J"'^ = IV(a35-H«5) f = Q'f (21) 

We see that the Eckart factorization is realized with as the charge density for the fluid. 
The equations of motion are given as 

D^,J''^' = 

FV = -iTT{a^g-'D^g) (22) 

where the color flow velocity u^' is defined by = nu^. We can also carry out a variation 
on g on the right, which leads to the equation d^j^ = 0. This equation is not independent of 
equations (22). The energy- momentum tensor for the fluid is given by 

T^^ = nF'ui'yy -g^'^inF'-F) (23) 
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and obeys the expected relation 

Q^rj^f^y ^ Tr( J^F^^^) (24) 

It is interesting at this stage to give an interpretation of the group element g. Let p = p"'t°' 
be the nonabelian charge density of a distribution of particles. Under a gauge transformation 
U G SU{2), it transforms as 

p' = U-^pU (25) 

We can diagonalize the hermitian matrix p at each point by an (x, i)-dependent transforma- 
tion. We may thus write p = g Pdiag 9~^j or 

= po Tria^g-H'^g) = n IV((735"'t^) (26) 

where pQ = n \s the eigenvalue of p. Evidently, n is gauge invariant. Comparing with our 
previous expression for = J"°, we see that the dynamical variable g{x, t) can be interpreted 
as the gauge transformation which diagonalizes the charge density at each point. The flow 
of the gauge invariant eigenvalues is given by u^. Under a gauge transformation, g — > U^^g. 

The generalization to higher groups is clear from this discussion. For a group G, we have 
rank{G) eigenvalues and hence rank{G) n's and u'^'s. The field g is an element of G and the 
action is given by 

rank{G) 

S = -i y ^''^ T^iQs g-^D^g) - J d^x F{ni, n2, . . .) + Sym (27) 

Qs are the diagonal generators of G, jg = n^. We have as many n's as there are simultaneously 
diagonalizable conserved charges. 

For a statistical distribution, we expect a chemical potential for each of the simultaneously 
diagonalizable conserved charges. The values of these chemical potentials are fixed by the 
values of the corresponding charges. Thus a statistical distribution of particles is specified 
by the values it has for the simultaneously diagonalizable conserved charges. What we find 
for the fluid is in conformity with this. 

One can give a general justification for (27) starting from a one-particle action again. The 
relevant observation is the following. The general form of KKS action is [10] 

S = -i^ws [ dt TtiQs g-' g) (28) 

There are quantization condition on this action. The numbers {wi, ■W2, ■ ■ ■) should form the 
weight vector of the highest weight state of a unitary irreducible representation of the group. 
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Upon quantization, the action (28) then gives a Hilbert space which corresponds to the 
unitary irreducible representation characterized by {wi,W2, ■ ■ ■)■ What we have given in (27) 

is indeed the appropriate generahzation of this result to fluids. (The action (28) leads to the 
symplcctic form which is the Kahler two- form on the coadjoint orbit of the element Wgqs- 
So this method is also known as the coadjoint orbit method.) 

The foregoing discussion also shows that ordinary hydrodynamics is a special case of this 
general structure where we have only one conserved charge, namely, the particle number. 

For the sake of completeness, we give the equations of motion for the action (27); these 

are 

df.j^ = 

Y,j^{D,Qsr = (29) 

s 

BF 

<^ = -iT^{qsg-^D^g) 

where = TY{(isg~^t"'g). The basic Poisson brackets are given by 

\p'^{x,tlp\y,t)] = r'>Yix,t)6^'\x-y) 

[jO(x,t),jO,(y,t)] = (30) 
[fsix,t),g{y,t)] = -ig{x) Qs 5^^\x - y) 

Here p" = is the color charge density and the color current is J"'^ = Ylsi'^Q's- 



3 Color skyrmions 



The description of nonabclian charge using the KKS form is a standard part of Lie group 
theory. In fact, if wc ask whether wc can find a classical theory which upon quantization 
gives finitc-dimcnsional Lie algebra (color) matrices in a single irreducible representation, 
the answer is the KKS action (28). Since the action for fluid dynamics which we postulate 
is a straightforward fluid generalization of this basic result for particles, we see that, quite 
generally, the color degrees of freedom of the fluid are described by g{x, t) e G and jg G 
Cartan subalgebra of G. For QCD, we have g{x,t) G SU{3). At a fixed time, we have 
g{x) : V — SU{3), where F is a region in in which the fluid exists. For configurations 
with 5 — 1 on the boundary of V, these functions g{x) are equivalent to g{x) : — SU{3). 
The homotopy classes, or equivalence classes up to smooth deformations, of such maps are 
given by n3[S'?7(3)] = Z. The topologically nontrivial configurations of g{x) are then solitons. 
These are mathematically the same as the usual skyrmions, although these are in the color 



8 



sector and not in the flavor sector [11, 12] ^. The soliton number which characterizes the 
homotopy classes n3[S'C/(3)], or the skyrmions, is 

Q = J d^x e^^^ T^{g-'d^g g-^d,g g-'d^g) (31) 

(This is essentially the color version of helicity.) 

For SU{3), there are two distinct types of maps g : ^ SU{3) wc can consider. One 
of them corresponds to the image of being an SU{2) subgroup of SU{3), the other corre- 
sponds to being mapped to the SO{3) subgroup defined by the generators (A2, — A5, A7), 
Aa being the usual Gell-Mann matrices of SU{3). Most of the features of the skyrmions we 
are considering will be clear from the first type of maps, namely, g : ^ SU{2) C SU{3), 
so we shall consider only these in this paper. This means that, effectively, we can restrict 
attention to SU{2) C SU{3). 

The explicit solution of the equation of motion to obtain the soliton will be very involved 
and will depend sensitively on the choice of Hamiltonian. But, as is usually done in the case 
of flavor skyrmions, we can choose an ansatz in a given topological sector, which depends 
on some parameters, and then variationally minimize the energy to fix the values of these 
parameters. In general, this will give a good qualitative (and, to some extent, even quanti- 
tative) description of the soliton. The simplest ansatz we can take for g{x) is the spherically 
symmetric ansatz 

gs{x) = cos (f){r) + i a ■ x sin^(r) (32) 

which leads to 

Q=^'°>-^'°°' (33) 

TT 

The profile of <^ as a function of r is not yet determined; for Q = 1, we need one step of 
height TT as we go from r = to r = 00. The simplest choice is the stereographic ansatz 

where R sets the scale for the soliton. In what follows we will make this simple choice. 

Up to this point we have not chosen a particular form for the Hamiltonian. The action 
(27) has the form 

S = -ij d^x f TT{a^g-^g) - J dt H (35) 

H is determined by the choice of F{n). We will take a simple form which is easy to work 
with. More specific choices, based on the equation of state of the color liquid, can be made. 



Skyrmions where the target space includes color have been considered in a different and unrelated context 
in [13]. 
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We do not expect the general features of the sohtons to be changed by this. Our choice for 
F{n) is then 

(36) 



J- J 



Ehminating j from the action (27), we identify 



2^2 



2 



(37) 



We see that, unhke the case of the flavor skyrmions, we will need an ansatz for as well. 
For this purpose, consider the collective coordinate quantization of the soliton gs{x). Color 
transformations act on the field g as g Ug. States with nonzero color charge can be 
generated from gs{x) by writing g{x,t) = U{t)gs{x). We also introduce spatial rotations 
R{t) G 50(3) and write a general ansatz as 



g{x,t) = U{t)gs{RS) 

= Uit)S-\t) g,{x) S{t) 



(38) 



where U{t), R{t) (or its (2 x 2)-matrix version S{t)) represent the collective coordinates. 
Notice that while the ansatz (32) is similar to what happens for flavor skyrmions, the collective 
coordinates enter into equation (38) in a very different way. Using (38) in the action, we find 



S = -iJdt[m Tr{a3U-^U) + (n - m) Tr((735-^5)] - J dt H 



(39) 



where m, n are given in terms of integrals involving combinations of j° and (f). The action 
is given in terms of the KKS forms for U and S and so lead to representations of the color 
algebra {SU{2) in this case) and the rotation algebra. For consistent quantization, with 
unitary representations of the two groups, we need m, n to be integers. Our ansatz for j° 
must be consistent with this. Since there are two conditions we take an ansatz for which 
has two functions in it; the simplest choice is 



(v ■ xY 



h{r) 



f = fir) + 

where Vi = iTr{a3g~^dig). Using the ansatz for g we then find 

j d^x f(r) = n 

J, 8hcj)' 



(40) 



/2 



15 



sm 



n — m 



The energy functional can be worked out to be 



S = 



2ir 



/a (, 



,9 2 sin"* 



+ 



2/x2 



d-^x 



(41) 



(42) 
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So far we have not used the specific form of (f). We now simphfy our ansatz further by taking 
h{r) to be a constant and /(r) = i(;(8(/>'^/15), where w is another constant. Further taking cf) 
to be given by the stereographic ansatz, the conditions (41) are evaluated as 



32tt'^R 



n 



15 



-w, 



n 



m 



15 



(w-h) 



which lead to 



15 



w = 



n. 



h = 



15 



327r2i? ' 327r2i? 
The energy function (42) can also be evaluated easily as 



(2m — n) 



£{R) 



S'k'^I?R + 



167rV2i?3 



{m — sf' + 



{m + s 



i2n 



(43) 



(44) 



(45) 



where s = jiP' ~ ^) ^Pi^ soliton and m defines the color charge. The collective 

coordinate action, upon quantization, gives an (m + l)-dimensional SU (2) color multiplet 
and (2s + l)-dimensional representation of the rotation group. The minimum of £{R) occurs 
at 

(m + s)2" 



{^iRT = 



15 1 
128^ 



(m - sY + 



and the variational estimate of the energy of the soliton is given as 



f = ^(30)V^ 



(m 



2^(rn + s)_ 



1/4 



(46) 



(47) 



Going back to (38), we see that the ansatz depends on R{t) G 50(3), so only the integer 
spins can arise. This means that the the difference n — m should be an even integer. With 
this condition, formula (47) gives the energy as a function of the color charge and spin. 



4 Discussion 

In this paper we considered the general formulation of nonabelian fiuid dynamics. Since it 
is based on symmetry considerations and the generalization of the mechanics of particles 
carrying nonabelian charges, there is a universality to this formulation. It is then a very 
general result that the theory admits color skyrmion configurations. We have obtained a 
(variational) mass formula for these as a function of spin and color charge, choosing a simple 
form for the Hamiltonian. Unlike the flavor skyrmions, there is no Wess-Zumino term, since 
there is no color anomaly. These solitons are bosonic excitations of integral spin. 

Clearly there are a number of interesting questions which need further investigation. 
There is a hierarchy of energy scales which is important in this context. The fluid approx- 
imation works in a regime where there are few single-particle hard scattering events, the 
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collective effects are important, and before rehadronization sets in. The form of the Hamil- 
tonian has to be fixed by the equation of state for the color degrees of freedom. (This is 
different from a general equation of state, since there are other conserved charges such as 
baryon number. Thus there are different partial pressures and here we need the equation of 
state for the color partial pressure as a function of the color charge density.) A more detailed 
investigation of how our results depend on the form of the Hamiltonian has to be carried 
out. Even after we have made the simple choice (37), the scale parameter is not specified. 
While it can be related to the equation of state, we expect that it is closely related to the 
magnetic screening in the plasma. One reason for this is to notice that, if we absorb the 
field 5 as a gauge transformation of A, in a way analogous to going to the unitary gauge in 
spontaneousy broken gauge theories, our Hamiltonian has the form H = J ^fi^A^ ■ + ■ ■ ■ . 

Another key issue is the production, detection and destruction of these solitons. The 
conservation of the topological charge is vitiated only if the field g{x, t) cannot be defined 
without singularities. So we expect that if such solitons are produced in nuclear collisions, it 
should happen at the transition to the deconfined fluid phase. Likewise, they should disappear 
due to rehadronization when the fluid description loses its validity. Dissipative effects can be 
introduced into our approach by modifying the equations of motion. While they can change 
many details, the qualitative features of the solitons should be unaltered. 

We hope to address some of the issues raised here in a future publication. 

VPN thanks I. Zahed for a useful comment. This work was supported in part by the National 
Science Foundation grant number PHY-0244873 and by a CUNY Collaborative Research 
Incentive grant. 
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